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Abstract. A new approach to the static route planning problem, based 
on a multi-staging concept and a scope notion, is presented. The main 
goal (besides implied efficiency of planning) of our approach is to 
address — with a solid theoretical foundation — the following two prac- 
tically motivated aspects: a route comfort and a very limited storage 
space of a small navigation device, which both do not seem to be among 
the chief objectives of many other studies. We show how our novel idea 
can tackle both these seemingly unrelated aspects at once, and may also 
contribute to other established route planning approaches with which 
ours can be naturally combined. We provide a theoretical proof that our 
approach efficiently computes exact optimal routes within this concept, 
as well as we demonstrate with experimental results on publicly available 
road networks of the US the good practical performance of the solution. 



1 Introduction 



The single pair shortest path SPSP problem in world road networks, also known 
as the route planning problem, has received considerable attention in the past 
decades. However classical algorithms like Dijkstra's or A* are in fact "optimal" 
in a theoretical sense, they are not well suitable for huge graphs representing 
real- world road networks having up to tens millions of edges. In such situation 
even an algorithm with a linear time complexity cannot be feasibly run on a 
mobile device lacking computational power and working memory. 

What can be done better? We focus on the static route planning problem 
where the road network is static in the sense that the underlying graph and 
its attributes do not change in time. Thus, a feasible solution lies in suitable 
preprocessing of the road network in order to improve both time and space 
complexity of subsequent queries (to find an optimal route from one position 
to another). However, to what extent such a preprocessing is limited in the size 
of the precomputed auxiliary data? It is not hard to see that there is always some 
trade-off between this storage space requirement and the efficiency of queries — 
obviously, one cannot precompute and store all the optimal routes in advance. 
See also a closer discussion below. 



Related Work. Classical techniques of the static route planning are rep- 
resented by Dijkstra's algorithm [pj, A* algorithm [IB] and their bidirectional 
variations [50]. In the last decade, two sorts of more advanced techniques have 
emerged and become popular. The first one prunes the search of Dijkstra's 
or A* algorithms using preprocessed information This includes, in particular, 
reach- based [15112] . landmarks |llll4j . combinations of those [T3], and recent 
hub-based labeling pQ. 

The second sort of techniques (where our approach conceptually fits, too) 
exploits a road network structure with levels of hierarchy to which a route can 
be decomposed into. For instance, highway and contraction hierarchies [21 24110] . 
transit nodes [3 , PCD [T5] and SHARC routing [4 1 5) represent this sort. Still, 
there are also many other techniques and combinations, but — due to lack of 
space — we just refer to Cherkassky et al. [5], Delling et al. [817] . and Schultes [2"3"] . 
Finally, we would like to mention the interesting notion of highway dimension 
[2], and the ideas of customizable [T] and mobile [22] route planning. 

Our Contribution. We summarize the essence of all our contribution already 
here, while we implicitly refer to the subsequent sections for precise definitions, 
algorithms, proofs, and further details. 

First of all, we mention yet another integral point of practical route plan- 
ning implementations — human-mind intuitiveness and comfortability of the com- 
puted route. This is a rather subjective requirement which is not easy to formalize 
via mathematical language, and hence perhaps not often studied together with 
the simple precise "shortest/fastest path" utility function in the papers. 

Intuitiveness and comfort of a route: Likely everyone using car navigation 
systems has already experienced a situation in which the computed route 
contained unsuitable segments, e.g. tricky shortcuts via low-category roads 
in an unfamiliar area. Though such a shortcut might save a few seconds on 
the route in theory, regular drivers would certainly not appreciate it and the 
real practical saving would be questionable. This should better be avoided. 

Nowadays, the full (usually commercially) available road network data con- 
tain plenty of additional metadata which allow it to detect such unreasonable 
routes. Hence many practical routing implementations likely contain some kinds 
of rather heuristic penalization schemes dealing with this comfortability issue. 
We offer here a mathematically sound and precise formal solution to the route 
comfort issue which builds on a new theoretical concept of scope (Sec. 

The core idea of a scope and of scope-admissible routes can be informally 
outlined as follows: The elements (edges) of a road network are spread into sev- 
eral scope levels, each associated with a scope value, such that an edge e assigned 
scope value s e is admissible on a route R if, before or after reaching e, such R 
travels distance less than s e on edges of higher scope level than e. Intuitively, 
the scope levels and values describe suitability and/or importance of particular 
edges for long-range routing. This is in some sense similar to the better known 
concept of reach [15] ; but in our case the importance of an edge is to be decided 



from available network metadata and hence its comfort and intuitive suitability 
is reflected, making a fundamental difference from the reach concept. 

The effect seen on scope admissible routes (Def. 13. f I I3.3l3.4p is that low-level 
roads are fine "near" the start or target positions (roughly until higher- level roads 
become available), while only roads of the highest scope levels are admissible in 
the long middle sections of distant routing queries. This nicely corresponds with 
human thinking of intuitive routes, where the driver is presumably familiar with 
neighborhoods of the start and the target of his route (or, such a place is locally 
accessible only via low- level roads anyway). On contrary, on a long drive the 
mentally demanding navigation through unknown rural roads or complicated 
city streets usually brings no overall benefit, even if it were faster in theory. 

To achieve good practical usability, too, road network segments are assigned 
scope levels (cf. Table [TJ according to expectantly available metadata of the net- 
work (such as road categories, but also road quality and other, e.g. statistical 
information). It is important that already experiments with publicly available 
TIGER/Line 2009 [19] US road data, which have metadata of questionable qual- 
ity, show that the restriction of admissible routes via scope has only a marginal 
statistical effect on shortest distances in the network. 

Furthermore, a welcome added benefit of our categorization of roads into 
scope levels and subsequent scope admissibility restriction is the following. 

Storage space efficiency: We suggest that simply allowing to store "linearly 
sized" precomputed auxiliary data (which is the case of many studies) may be 
too much in practice. Imagine that setting just a few attributes of a utility 
function measuring route optimality results in an exponential number of 
possibilities to which the device has to keep separate sets of preprocessed 
auxiliary data. In such a case a stricter storage limits should be imposed. 

In our approach, preprocessed data for routing queries have to deal only with 
the elements of the highest scope(s). This allows us to greatly reduce the amount 
of auxiliary precomputed information needed to answer queries quickly. We use 
(Sec. 14. 1|) a suitably adjusted fast vertex-separator approach which stores only 
those precomputed boundary-pair distances (in the cells) that are admissible 
on the highest scope level. This way we can shrink the auxiliary data size to 
less than 1% of the road network size (Table [3]) which is a huge improvement. 
Recall that vertex-separator preprocessing produces a partition of the network 
graph into moderate-size cells (several thousands of edges, say) such that only 
(selected) distances between pairs of their boundary vertices are precomputed. 

Not to forget, our subsequent routing query algorithm (Sec. [5]) then answers 
quickly and exactly (not heuristically) an optimal route among all scope admis- 
sible ones between the given positions. The way we cope with scope admissibility 
(among other aspects) in a route planning query, using the precomputed auxil- 
iary data, It is briefly summarized as follows. 

Multi-staging approach: The computation of an optimal route is split 
into two (or possibly more - with finer network metadata and hierarchi- 
cal separators) different stages. In the local - cellular, stage, a modification 



of plain Dijkstra's algorithm is used to reach the cell boundaries in such a 
way that lower levels are no longer admissible. Then in the global - bound- 
ary, stage, an optimal connection between the previously reached boundary 
vertices found on the (much smaller) boundary graph given by auxiliary 
data. 

Notice that the cellular stage may possibly cross the boundaries of a few 
adjacent cells if the start or target is near to them, but practical experiments 
show that such a case is quite rare. After all, the domain of a cellular stage is 
a small local neighbourhood, and Dijkstra's search can thus be very fast on it 
with additional help of a reach- like parametrization (Def. !3-8[) . Then, handling 
the precomputed boundary graph in the global stage is very flexible — since no 
side restrictions exist there — and can be combined with virtually any other es- 
tablished route planning algorithm (see Sec. [5]). The important advantage is that 
the boundary graph is now much smaller (recall, < 1% in experiments) than the 
original network size, and hence computing on it is not only faster, but also more 
working-memory efficient which counts as well for mobile navigation devices. 

Paper organization. After the informal outline of new contributions, this pa- 
per continues with the relevant formal definitions — Section [5] for route planning 
basics, and Section [3] for thorough description of the new scope admissibility 
concept. An adaptation of Dijkstra for scope is sketched in Sec. 13.11 Then Sec- 
tion @] shows further details of the road network preprocessing (|4.1[) and query 
(O algorithms. A summary of their experimental results can be found in Tables 
[3] and [6] 

2 Preliminaries 

A directed graph G is a pair of a finite set V(G) of vertices and a finite multiset 
E{G) C V{G) x V(G) of edges. The reverse graph G R of G is a graph on the 
same set of vertices with all of the edges reversed. A subgraph H of a graph 
G is denoted by H C G. A subgraph H C G is induced by a set of edges 
F C E(G) if E(H) = F and V(H) = {u e V(G) | 3/ 6 F incident with u}; we 
then write H = G[F]. 

A walk P C G is an alternating sequence (uq, e±, Mi, . . . , e k , Uk) of vertices and 
edges of G such that = (iti_i, Ui) for i = 1, . . . , k. A subwalk is a subsequence 
of a walk. A concatenation P\. P2 of walks P\ = (uq, &i, u±, . . . , e kl Uk) and P2 = 
(u fc ,e fc+ i,u fc+ i, . . . ,ei,ui) is the walk (u ,ei,Mi, . . . , e k ,u k , e k +\, . ■ ■ , e h u x ). If P 2 
is a single edge /, then we write Pi. f. 

A walk Q is is a prefix of another walk P if Q is a subwalk of P start- 
ing with the same index, and analogically with suffix. The prefix set of a walk 
P = (u ,e 1 ,...,e k ,u k ) is Prefix(P) = {(u , e x , . . . , e h Ui)\ < i < fc}, and 
analogically Suffix(P) = {(u.i, e^+i, . . . , e k , u k )\ < i < fc}. Two walks are over- 
hanging (one another) if either one is a subwalk of the other, or a non-zero-length 
suffix of one is a prefix of the other (informally, one can traverse both with one 
superwalk) . 



The weight of a walk P C G wrt. a weighting w : E(G) i— J> M of G is defined 
as |P| M , = w(ei) + t«(e2) + - • •+w(ei e ) where P = (uq, e±, . . . , et,Uk)- The distance 
5 w (u, v) from it to t> in G is the minimum weight over all walks from u to v, or 
oo if there is no such walk. 

A road network can be naturally represented by a graph G such that the 
junctions are represented by V(G) and the roads (or road segments) by E{G). 
Chosen cost function (e.g. travel time, distance, expenses, etc.) is represented by 
a non-negative edge weighting w : E(G) >-> R^" of G. 

Definition 2.1 (Road Network). Let G be a graph with a non-negative edge 
weighting w. A road network is the pair (G,w). 

A brief overview of classical Dijkstra's and A* algorithms and their bidirec- 
tional variants for shortest paths follows, but we also would like to recall the 
useful notion of a reach given by Gutman |15j . 

Definition 2.2 (Reach [15j). Consider a walk P in a road network (G,w) 
from s to t where s, t 6 V(G). The reach of a vertex v 6 V(P) on P is rp(v) = 
min{\P sv \ w , \P vt \ w } where P sv and P vt is a subwalk of P from s to v and from 
v to t, respectively. The reach of v in G, r(v), is the maximum value of tq(v) 
over all optimal walks Q between pairs of vertices in G such that v € V{Q). 

Classical Shortest Paths Algorithms. Classical Dijkstra's algorithm solves 
the single source shortest paths problem^ in a graph G with a non-negative 
weighting w. Let s € V(G) be the start vertex (and, optionally, let t 6 V(G) be 
the target vertex). 

— The algorithm maintains, for all v € V{G), a (temporary) distance estimate 
of the shortest path from s to v found so far in d[v], and a predecessor of v 
on that path in ir [v] . 

— The scanned vertices, i.e. those with d[v] = 5 w (s, v) confirmed, are stored in 
the set T; and the discovered but not yet scanned vertices, i.e. those with 
oo > d[v] > S w (s,v), are stored in the set Q. 

— The algorithm work as follows: it iteratively picks a vertex u £ Q with 
minimum value d[u] and relaxes all the edges (u, v) leaving u. Then u is 
removed from Q and added to T. Relaxing an edge (u, v) means to check if 
a shortest path estimate from s to v may be improved via u; if so, then d[v] 
and tt[v] are updated. Finally, v is added into Q if is not there already. 

— The algorithm terminates when Q is empty (or if t is scanned). 

Time complexity depends on the implementation of Q\ such as it is 
0{\E(G)\ + \V{G)\ log |V(G)|) with the Fibonacci heap. 

A* algorithm is also known as goal directed search and it is equivalent to 
aforementioned Dijkstra's one on a graph modified as follows: 

1 Given a graph and a start vertex find the shortest paths from it to the other vertices. 



— A potential function p : V(G) <— > R is defined and used to reduce edge 
weights such that w p (u, v) — w(u, v) —p(u) +p(v). If w p is non-negative, the 
algorithm is correct and p is called feasible. 

— Each path from s to t found by Dijkstra's algorithm in G with the reduced 
weighting w p then differs from reality by pit) — p(s) and this value must be 
subtracted at the end. 

Dijkstra's and A* algorithms can be used "bidirectionally" to solve SPSP 
problem. Informally, one (forward) algorithm is executed from the start vertex 
in the original graph and another (reverse) algorithm is executed from the target 
in the reversed graph. Forward and reverse algorithms can alternate in any way 
and algorithm terminates, for instance, when there is a vertex scanned in both 
directions. 

Dijkstra's and A* algorithms can be accelerated by reach as follows: when dis- 
covering a vertex v from u, the algorithm first tests whether r(v) > d[u] + w(u, v) 
(the current distance estimate from the start) or r(v) > loweriv) (an auxiliary 
lower bound on the distance from v to target), and only in case of success it 
inserts v into the queue of vertices for processing. 

3 The New Concept — Scope 

The main purpose of this section is to provide a theoretical foundation for the 
aforementioned vague objective of "comfort of a route" . Recall that the scope 
levels referred in Definition 13.11 are generally assigned according to auxiliary 
metadata of the road network, e.g. the road categories and additional available 
information which is presumably included with it; see Table [T] Such a scope level 
assignment procedure is not the subject of the theoretical foundation. 

Definition 3.1 (Scope). Let (G,w) be a road network. A scope mapping is 
defined as S : E{G) H> No U {oo} such that 0, oo G Im(S). Elements of the 
image Im{S) are called scope levels. Each scope level i G Im(S) is assigned a 
constant value of scope vf G KoU{oo} such that = Vq < vf < • • • < = oo. 

To give readers a better feeling of how the scope level assignment outlined in 
Table [T] looks like, we present some statistics of the numbers of edges assigned 
to each level in Table [5] 

There is one more formal ingredient missing to make the scope concept a 
perfect fit: imagine that one prefers to drive a major highway, then she should 
better not miss available slip-roads to it. This is expressed with a "handicap" 
assigned to the situations in which a turn to a next road of higher scope level is 
possible, as follows: 

Definition 3.2 (Turn-Scope Handicap). Let S be a scope mapping in (G, w). 
The turn-scope handicap hg{e) G No U {oo} is defined, for every e = (u,v) G 
E(G), as the maximum among S(e) and all S(f) over f — (v,w) G E(G). Each 
handicap level i is assigned a constant nf such that < Kq < ■ ■ ■ < . 



Table 1. A very simple demonstration of a scope mapping which is based just on 
the road categories. Highways and other important roads have unbounded (oo) scope, 
while local, private or restricted roads have smaller scope. The zero scope is reserved 
for roads that physically cannot be driven through (including, for instance, long-term 
road obstructions). The weight function of this example is the travel distance in meters. 



Scope level i 


? — 

Value vf 


o — 

Handicap K i 


Road category 








1 


Alley, Walkway, Bike Path, Bridle Path 


1 


250 


50 


Parking Lot Road, Restricted Road 


2 


2000 


250 


Local Neighborhood Road, Urban Roads 


3 


5000 


600 


Rural Area Roads, Side Roads 


oo 


oo 


(oo) 


Highway, Primary (Secondary) Road 



The desired effect of admitting low-level roads only "near" the start or target 
positions — until higher level roads become widely available — is formalized in next 
Def. 13.31 13.41 We remark beforehand that the seemingly complicated formulation 
is actually the right simple one for a smooth integration into Dijkstra. 

Definition 3.3 (Scope Admissibility). Let (G,w) be a road network and let 
x G V(G). An edge e — (u, v) G E(G) is X- admissible in G for a scope mapping 
S if, and only if, there exists a walk P C G — e from x to u such that 

1. each edge of P is x-admissible in G — e for S, 

2. P is optimal subject to (1), and 

3. for £ = S(e), J2feE(P),s(j)>e w (f) + ^feE(P),h s (f)>e>s(f) K f — v t ■ 

Note; every edge e such that 5(e) = oo (unbounded scope level) is always 
admissible, and with the values of vf growing to infinity, Def. l3.3l tends to admit 
more and more edges (of smaller scope). 

Definition 3.4 (Admissible Walks). Let (G,w) be a road network and S a 
scope mapping. For a walk P = (s = uq, e±, . . . e/., Uk — t) CG from s to t; 

— P is s-admissible in G for S if every ej £ E(P) is s-admissible in G for S , 

— and P is st-admissible in G for S if there exists < j < k such that every 
&i G E(P), i < j is s-admissible in G, and the reverse of every G E(P), 
i > j is t-admissible in G R - the graph obtained by reversing all edges. 



Table 2. An example of the scope levels distribution in several road networks graphs 
of the US; the initial scope is assigned in accordance with road categories as described 
in Table [T] and then balanced algorithmically. 



Scope level 


US A- all 


Alabama 


Connecticut 


Florida 


Georgia 


Indiana 


oo 


8.171% 


9.726% 


12.874% 


8.392% 


11.331% 


8.672% 


3 


62.420% 


73.711% 


64.022% 


58.389% 


69.564% 


64.967% 


2 


19.712% 


10.052% 


17.018% 


28.687% 


13.250% 


17.270% 


1 


6.570% 


5.933% 


3.978% 


2.038% 


4.946% 


7.239% 





0.660% 


0.018% 


0.092% 


0.083% 


0.020% 


0.007% 



Proper Scope Mapping. In a standard connectivity setting, a graph (road 
network) G is routing- connected if, for every pair of edges e, f £ E(G), there 
exists a walk in G starting with e and ending with /. This obviously important 
property can naturally be extended to our scope concept as follows. 

Definition 3.5 (Proper Scope). A scope mapping S of a routing- connected 
graph G is proper if, for all i g Im(S), the subgraph G^ induced by those edges 
e G E(G) such that 5(e) > i is routing- connected, too. 

Note that validity of Definition 13.51 should be enforced in the scope- 
assignment phase of preprocessing (e.g., the assignment should reflect known 
long-term detours on a highwajQ accordingly) . The topic of connectivity require- 
ments for assigned scope levels deserves a bit closer explanation. A good scope 
mapping primarily reflects the road quality (given by its category and other 
attributes) as given in the road network metadata. This simple assignment is, 
however, not directly usable in practice due to rather common errors in road 
network data. Such errors typically display themselves as violations of connec- 
tivity at different scope levels, i.e. as violations to Def. 13.51 Hence these errors 
can be algorithmically detected and subsequently repaired (preferably with help 
of other road network metadata) so that the final scope mapping conforms to 
Def. 13.51 -being a proper scope. 

Existence of admissible walks. A natural conclusion of routing connectivity 
(proper scope) of a road network then reads 

Theorem 3.6. Let (G, u>) be a routing- connected road network and let S be 
a proper scope mapping of it. Then, for every two edges e — (s,x) G E{G), 
f = (y,t) G E(G), there exists a st-admissible walk P C G for S such that P 
starts with the edge e and ends with f. 

Proof. We refer to the notation of Definition 13.51 Let i — min(6>(e), S(f)), and 
so e, / G . Since S is proper, there exists a walk P l C GN starting with the 
edge e and ending with /, and \P l \ w < oo. If i = oo, then we are done since every 
edge of G*-°°' is automatically st-admissible (regardless of s,t). So, by means of 
contradiction, we assume that i e Nfl Im(S) is highest possible such that the 
theorem failed, i.e. no such P % is st-admissible. 

Let P l — (i>o,ei, . . . ,ek,Vk). We find maximal index 1 < p < k such that 
<5( e p) > « + 1 or that hs(e p ) > i + 1 (where the value of hs(e p ) is witnessed by 
an edge /o of <S(/o) > « + 1 starting in v p ). Informally, p is the least position 
on P l where the scope admissibility condition gets affected at level i — see the 
formula in Def. 13.31 Since P % is not st-admissible by our assumption, this p is 
well defined. In the first case, S(e p ) > i + 1, we set e' = e p . In the second case, 
S(e p ) < i and <S(/o) > i + 1, we set e' — fo. We also denote by s' the tail of e'. 
We symmetrically define /' and t' in the reverse walk of P % in G R . 

2 Note, regarding a real-world navigation with unexpected road closures, that the pro- 
per-scope issue is not at all a problem — a detour route could be computed from the 
spot with "refreshed" scope admissibility constrains. Here we solve the static case. 



Now j — min(<S(e'), <S(/')) > Since S is proper, there exists a walk P 3 C 
G^l starting with the edge e' and ending with /', and \P^\ W < oo. Moreover, 
by maximality of our choice of i, this P- 3 can be chosen such that P J is s't'- 
admissible. Then we define P' C P % U P J as the walk which starts and ends as 
P l while using P- 7 in the section "between" e' and /'. By the least choice of 
e' (/') above, all edges of P' are clearly also si-admissible (which follows from 
admissibility of all edges of P J at levels > i there). The proof is finished. □ 

Well distributed scope. There is another, more vague, technical requirement 
on a useful scope mapping, which becomes particularly important in the opening 
cellular stage of the query algorithm in Section [5] for each bounded scope level 
there should be no relatively large areas (subnetworks) not penetrated by any 
road of higher scope level. This requirement is formally described as that the 
subgraph G — V(G™), i.e. the subgraph induced by those vertices incident only 
with edges of scope level < i, has no relatively large connected components for 
each i g Im(S). Such as, for i — oo the words "no relatively large" mean no 
component of size significantly larger than the expected cell size. For smaller 
values of i the size limit on components is accordingly smaller. We then say that 
the scope mapping S is well-distributed in G. 

Again, as in the case of proper scope, the requirement of having well- 
distributed scope mapping is not strictly a subject of the theoretical founda- 
tion of scope. Instead, wisely designed road networks (with the corresponding 
metadata) should conform to such a requirement automatically. In other words, 
if a scope level assignment does not produce a well-distributed scope mapping, 
then there is something wrong; either directly in the network design, or in the 
provided metadata (or in the way we understand it). Yet, violations of the well- 
distributivity property can be easily detected in an algorithmic way, and also au- 
tomatically repaired via taking some heuristically selected shortest walks across 
the large component, and "upgrading" them into higher scope levels. 

3.1 <S-Dijkstra's Algorithm and 5-Reach 

As noted beforehand, Def. 13.41 smoothly integrates into bidirectional Dijkstra's 
or A* algorithm, simply by keeping track of the admissibility condition (3.): 

iS-Dijkstra's Algorithm (one direction of the search). 

— For every accessed vertex v and each scope level I € Im(S), the algorithm 
keeps, as <Je[v], the best achieved value of the sum. formula in Def. 13.31 (3.). 

— The s-admissibility of edges e starting in v depends then simply on o\s( e ) [ v ] < 
f§(e) ' an d only s-admissible edges are relaxed further. 

The full details follow here, in Algorithm 13.11 Notice that, although a route 
planning position in a road network is generally an edge (not just a vertex) 
there is no loss of generality if we consider only vertices as the start and target 
positions of a route planning query. In the general case when a position is an 
edge /, we may simply subdivide / 



Algorithm 3.1 5-Reach-based 5-Dijkstra's Algorithm 



Input: A road network (G,w), a scope S, a start vertex s G V(G), an 5-reach r s . 
Output: For every v G V(G), an optimal s-admissible walk from s to v in G (or oo). 
Relax(u, v, 7) 

// Temporary distance estimate updated. 



2 
3 
4 
5 
6 
7 
8 
9 
10 



if d[u] + w(u, v) < d[v] then 
Q^Qu{v} 

d[v] <— d[u] +w(u,v); tt[v] <— u 

fi 

if d[u] + w(u,v) < d[v] then 
for all i G Im(S) do 

<7ik>] min{cri[«], o~% [u] + 7i} 
done 

fl 

return 



// Scope admissibility vector updated. 



5-Dijkstra(G, w, 5, s, r 5 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 



for all v G V(G) do 

d[v] <— 00; tv[v] <— _L; 

a[u] (00, . . . , 00) 
done 

d[s] <- 0; Q <- {s}; cr[s] <- (0, . . . , 0) 

while Q / do 

w min d [](Q); Q «- Q \ {u} 
for all / = (it, v) G E(G) do 



// Initialization. 
/ / Distance estimate and predecessor. 

// Scope admissibility vector. 



// Main loop processing all vertices. 
// Pick a vertex u with the minimum d[u]. 

// All edges from u; subject to 



if r*(v) > d[u]+w(f) V S(f) = 00 then 



if 



// S -reach check, and 
°~s(f)[ u ] 5: u s(f) then // s- admissibility check. 

for all i G Im{S) do // Adjustment to scope admissibility. 
if S(f) > i then 7, w(f) else 7^ ^— fi 



if h s (f) > i > «S(/) then 7l 
done 

Relax(u, 11,7) 



<- 7i + «f fi 

// Relaxation of f = (it, «) 



fi 

done 
done 

ConstructWalk (G, d, 7r) 



// Postprocessing - generating output. 



Theorem 3.7. 5-Dijkstra's algorithm, for a road network (G, w), a scope map- 
ping S, and a start vertex s G V(G), computes an optimal s-admissible walk 
from s to every v G V(G) in time 0{\E(G)\ ■ \Im(S)\ + \V(G)\ ■ \og\V(G)\). 

Proof. We divide the proof of correctness of the algorithm into three steps: 

1. Assume that — in line 10 — the algorithm correctly identifies the s-admissible 
edges /, and that — in line 9 — no edge / is discarded which is a part of some 
optimal s-admissible walk between the start s and any other s-saturated 



vertex y (in G for S). Then we can directly use a traditional proof of or- 
dinary Dijkstra's algorithm to argue that the computed results are optimal 
s-admissible walks starting from s in the network (G, w) for scope S. 
As to the claimed complexity bound, now it is enough to add a factor of 
|7m(iS)| (which can be regarded as a constant) for the loop on line 11. 

2. Suppose the assumption on line 9 got potentially violated, i.e. the particular 
edge / = (u, v) of bounded scope belongs to an optimal s-admissible walk Q 
from s to some s-saturated y. Let Q' C Q be a subwalk starting with x — s 
and ending with the edge / (i.e., in the vertex v). By Definition 13.81 the 
S-reach of v is defined and at least \Q'\ W > d[u] +w(f), and so the condition 
on line 9 is evaluated as TRUE. The claim is done. 

3. Hence it is enough to prove that the scope- admissibility vector a is computed 
such that an edge / = (u, v) starting in a vertex u is s-admissible if and 
only if Gs{j)\u\ < v s(f)' ^ Definition 13.31 and induction on the length 
(number of edges) of the discovering walk from s to u; it is just enough to 
straightforwardly verify that the adjustments to er accumulated in lines 12- 
13 exactly correspond to the summation formula in Definition 13.31 (3). and 
employ the known fact that d[u] is an optimal weight of a walk to u. □ 

Furthermore, practical complexity of this algorithm can be largely decreased 
by a suitable adaptation of the reach concept (Def. I2.2[) . given in Def. 13.81 For 
x G V(G) in a road network with scope S, we say that a vertex u G V(G) is 
x- saturated if no edge / = (u, v) of G from u of bounded scope (i.e., S(f) < oo) 
is x-admissible for S. A walk P with ends s, t is saturated for S if some vertex 
of P is both s-saturated in G and t-saturated in the reverse network G R . 

Definition 3.8 (5-reach). Let (G,w) be a road network and S its scope map- 
ping. The 5-reach of v G V(G) in G, den. r s (v), is the maximum value among 
\P XV \ W over all x, y G V(G) such that y is x-saturated while v is not x-saturated, 
and there exists an optimal x- admissible walk PCG from x to y such that P xv 
is a subwalk of P from x to v. r (v) is undefined (oo) if there is no such walk. 

There is no general easy relation between classical reach and 5-reach; they 
both just share the same conceptual idea. Moreover, iS-reach can be computed 
more efficiently (unlike reach) since the set of non- x-saturated vertices is rather 
small and local in practice, and only its x-saturated neighbors are to be consid- 
ered among the values of y in Def. 13.81 

The way 5-reach of Def. 13.81 is used to amend 5-Dijkstra's algorithm is again 
rather intuitive; an edge / = (u, v) G E(G) is relaxed from u only if S(f) = oo, 
or r s (v) > d[u] + w(f). 

Theorem 3.9. Assume 5-Reach-based 5-Dijkstra's algorithm (Ala. \3.1}) . with 
a road network (G,w), a scope mapping S, and any upper bound r s on the S- 
reach in G. This algorithm, for a start vertex s G V(G) as its input, computes, 
for every v G V(G), an optimal s-admissible walk from s to v in G in time 
0{\E{G)\-\Im{S)\ + \V{G)\.\og\V{G)\). 



Remark 3.10. Notice that we have used an "upper bound on the iS-reach" in 
the statement, although we can easily compute the 5-reach exactly. The deeper 
reason for doing this generalization is that in a practical case of multiple utility 
weight functions for the same road network, we may want to store just one 
maximal instance among the iS-reach values for every vertex in the network 
(instead of a separate value per each utility function). 

However, our 5-reach amending scheme has one inevitable limit of usability — 
it becomes valid only if the both directions of Dijkstra's search get to the "sat- 
urated" state. (In the opposite case, the start and target are close to each other 
in a local neighborhood, and the shortest route is quickly found without use of 
iS-reach, anyway.) Hence we conclude: 

Theorem 3.11. Let s,t 6 V(G) be vertices in a road network (G,w) with a 
scope mapping S. Bidirectional 5-reach 5-Dijkstra's algorithm computes an op- 
timal one among all st-admissible walks in G from s to t which are saturated 
for S. 

Notice that this strengthening Theorem 13.111 is directly implied by Theo- 
rem 133] which also directly implies former Theorem 13 .71 fwhen setting r s = oo). 

4 The Route-Planning Algorithm — Preprocessing 

Following the informal outline from the introduction, we now present the second 
major ingredient for our approach; a separator based partitioning of the road 
network graph with respect to a given scope mapping. 

4.1 Partitioning into Cells 

At first, a road network is partitioned into a set of pairwise edge-disjoint sub- 
graphs called cells such that their boundaries (i.e., the vertex-separators shared 
between a cell and the rest) contain as few as possible vertices incident with 
edges of unbounded scope. The associated formal definition follows. 

Definition 4.1 (Partitioning and Cells). Let£ = {E\, . . . ,Eg} be a partition 
of the edge set E(G) of a graph G. We call cells of {G,£) the subgraphs O; = 
G[Ei] C G, for i = 1,2, ... , i. The cell boundary r(Cj) of C t is the set of all 
vertices that are incident both with some edge in Ei and some in E(G) \ E i} and 
the boundary of £ is r(G,£) = \Ji<i<t r ( C i)- 

Practically, we use a graph partitioning algorithm hierarchically computing 
a so-called partitioned branch- decomposition of the road network. The algorithm 
employs an approach based on max-flow min-cut which, though being heuristic, 
performs incredibly well — being fast in finding really good small vertex separa- 
torsJl 

3 It is worth to mention that max-flow based heuristics for a branch-decomposition 
have been used also in other combinatorial areas recently, e.g. in the works of Hicks. 



Partitioning Algorithm. Notice that we are decomposing the whole graph G 
and not only its unbounded-scope subgraph G^ 00 ' , this is because of the cellular 
stage of our query algorithm. A brief outline of the partitioning method follows. 

1. Simplification step. At first, all weakly disconnected components are placed 
into a special "disconnected" cell and then removed from the road network to 
ensure that the rest is weakly-connected. Then the road network is contracted 
as much as possible so that all self-loops and parallel edges are removed, 
each maximal induced subtree is cut away, all maximal induced subpaths 
are contracted into single edges. This way we reduce the size of the road 
network by 10%, approximately. 

2. Replacement step. We are interested in vertex cuts and thus, prior to max- 
flow min-cut computation, we have to replace vertices by new edges as fol- 
lows: each vertex v in the road network is replaced by a new edge (v{, e v ,v ) 
such that every edge with its end in v has now its end in u, and every edge 
with its start in v has now its start in v Q . Notice that any cut consisting of 
such new edges represents a vertex cut in the original road network. 

3. Capacity assignment. In general, each newly introduced edge has unit ca- 
pacity and capacities of the other edges are set to positive infinity so that 
they will never be saturated during the max-flow computation and hence 
only new edges could appear in a cut (so that they form a vertex cut in the 
original road network). It is also beneficial (to allow a certain level of par- 
allelism) to cut the graph into two subgraphs of similar sizes and hence we 
choose sufficiently distant vertices to be the source and the sink and, more- 
over, every edge within some "restricted distance" from the source or sink 
gets infinite capacity to prevent its saturation (this way we can ensure that a 
cut appears somewhere in the "middle" between the source and sink). Also, 
to minimize the size (the number of edges) of the strict boundary graph we 
add some handicap to the capacities of those new edges which are adjacent 
to network edges of unbounded scope level. 

4. Max-flow min-cut computation. As soon as all capacities are assigned, the 
maximum flow is computed and we obtain a minimum vertex cut in the 
original road network. Hence we cut the original road network into two sub- 
networks (of similar size), i.e. into two super-cells. The whole process now 
continues from step (1) concurrently in these subnetworks until all cells have 
suitable size. 

Furthermore, our algorithm always leverages natural disposition of a road 
network such that each sufficiently big autonomous region (i.e. the US state) is 
partitioned separately. The resulted partitions are then "glued" together using 
natural borders. Clearly, this process might create too big boundaries and then 
the algorithm unions adjacent borderline cells and tries to partition them bet- 
ter. Such a separate partitioning makes it possible to employ a very beneficial 
parallelism. 



4.2 Bounday Graph Construction 

Secondly, the in-cell distances between pairs of boundary vertices are precom- 
puted such that only the edges of unbounded scope are used. This simplification 
is, on one hand, good enough for computing optimal routes on a "global level" 
(i.e., as saturated for scope in the sense of Sec. I3.1[) . On the other hand, such a 
simplified precomputed distance graph (cf. Bg) is way much smaller than if all 
boundary-pair distance were stored for each cell. See in Tabled! the last column. 

We again give the associated formal definition and a basic statement whose 
proof is trivial from the definition. 

Definition 4.2 (Boundary Graph; 5-restricted). Assume a road network 
(G, w) together with a partition £ = {E\, . . . , Eg} and the notation of Def. 
For a scope mapping S of G, let G^ 00 ' denote the subgraph of G induced by the 
edges of unbounded scope level oo, and let G°° = G[Ej\ fl G^. 

The (5-restricted) in-cell distance graph Di of the cell Ci is defined on the 
vertex set V{D{) = r{Ci) with edges and weighting pi as follows. For u, v £ 
r(Ci) n V(C°°) only, let 8^{u,v) be the distance in G?° from u to v, and let 
f = (u,v) e E(Di) iffpi(f) := 6%{u,v) < oo. 

The weighted (5-restricted) boundary graph Bg of a road network (G, w) 
wrt. scope mapping S and partition £ is then obtained as the union of all the 
cell-distance graphs Di for i = 1, 2, . . . simplified such that for each bunch of 
parallel edges only one of the smallest weight is kept in Bg . 

Proposition 4.3. Let (G, w) be a road network, S a scope mapping of it, and 
£ a partition of E(G). For any s,t £ r(G,£), the minimum weight of a walk 
from s to t in G^°°^ equals the distance from s to t in Bg . 

4.3 Experimental Evaluation. 

The prototype of the preprocessing algorithm is written in C and uses Ford- 
Fulkerson's max-flow min-cut algorithm and cells of approximately 5000 edges. 
Minimum possible cell size is hard-coded to 2000 and maximum to 10000 edges. 
Publicly available road networks are taken from TIGER/Line 2010 [TH] pub- 
lished by the U.S. Census Bureau using directed edges, i.e. every traffic lane is 
represented by one edge. The compilation is done by gcc 4.3.2 with -02, and 
the preprocessing has been executed on a quad-core XEON machine, with 16GB 
RAM running Debian 5.0.4, GNU/Linux 2.6.26-2-xen-amd64. 

We have run our partitioning algorithm (Def. 14. ip , in-cell distance computa- 
tions (Def. I4.2[) . and 5-reach computation (Def. I3.8[) in parallel on a quad-core 
XEON machine with 16 GB in 32 threads. A decomposition of the continental 
US road network into the boundary graph Bg, together with computations of 
5-reach in G (and iS-reach in G R ), distances in DiS, and standard reach estimate 
for Bg, took only 192 minutes altogether. This, and the tiny size of Bg, are both 
very for potential practical applications in which the preprocessing may have 
to be run and the small boundary graphs separately stored for multiple utility 
weight functions (while the S- reach values could still be kept in one maximizing). 



The collected data are briefly summarized in Table 0] and Table [3J To get a 
taste of the topic, some examples of local cell partitioning results are depicted 
in Fig. H and [TJ 



Table 3. Partitioning results for the TIGER/Line 2009 [19] US road network. The left 
section identifies the (sub)network and its size, the middle one the numbers and average 
size of partitioned cells, and the right section summarizes the results — the boundary 
graph size data, with percentage of the original network size. This boundary graph 
Be is wrt. the simple scope assignment of Table |3~T1 and V(Bs) includes also isolated 
boundary vertices (i.e. those with no incident edge of unbounded scope) . Notice the tiny 
size (< 1%) of Be compared to the original road network, and the statistics regarding 
cell boundary sizes: maximum is 74, average 19, median 18, and 9-decil is 31 vertices. 



Input 


G 








Partitioning 


Boun 


dary Graph 


B £ 


Road network 


#Edg 


;es 


#Cells 


Cell sz. 


^Vertices 


#Ed 


ges / % size 


USA-all 


88 


742 


994 


15 


862 


5 594 


253 641 


524 


872 


1 0.59% 


USA-east 


24 


130 


457 


4 


538 


5 317 


62 692 


107 


986 


/ 0.45% 


USA-west 


12 


277 


232 


2 


205 


5 567 


23 449 


42 


204 


/ 0.34% 


Texas 


7 


270 


602 


1 


346 


5 366 


17 632 


36 


086 


/ 0.50% 


California 


5 


503 


968 


1 


011 


5 444 


11 408 


16 


978 


/ 0.31% 


Florida 


3 


557 


336 




662 


5 373 


5 599 


25 


898 


/ 0.73% 




Fig. 2. The partition computed using a max- flow algorithm respects the natural road 
network disposition so that cells boundaries often lie along the rivers, railways, canyons, 
etc. This, in particular, leads to smaller boundaries. New Haven, CT, United States. 



Fig. 1. The road network partitioned into two cells (blue and green) with boundaries 
of size two (red), notice that each walk starting in a cell and ending in another must 
pass the boundary. Acadia National Park, Mount Desert, ME, United States. 



Table 4. The example of the partitioning of several US states. 



Road network 


# Vertices 


#Edges 


#Cells 


Aug. cell size 


Time 


Alabama 


811 434 


1 


926 


052 


338 


5 


698 


11 


min 


Indiana 


800 295 


1 


975 


898 


372 


5 


311 


16 


min 


Michigan 


860 421 


2 


145 


960 


399 


5 


378 


61 


min 


Minnesota 


908 292 


2 


166 


138 


383 


5 


666 


18 


min 


Arizona 


893 163 


2 


184 


866 


420 


5 


202 


25 


min 


Georgia 


945 212 


2 


226 


392 


400 


5 


655 


14 


min 


New York 


890 684 


2 


236 


530 


422 


5 


299 


47 


min 


North Carolina 


1 104 258 


2 


497 


764 


409 


6 


107 


12 


min 


Oklahoma 


1 049 680 


2 


508 


862 


436 


5 


754 


23 


min 


Ohio 


1 085 287 


2 


761 


920 


483 


5 


404 


38 


min 


Illinois 


1 085 287 


2 


761 


920 


530 


5 


211 


21 


min 


Missouri 


1 291 751 


3 


020 


152 


519 


5 


891 


34 


min 


Pennsylvania 


1 263 737 


3 


081 


096 


574 


5 


367 


26 


min 


Virginia 


1 489 661 


3 


333 


864 


479 


6 


960 


20 


min 


Florida 


1 387 005 


3 


488 


194 


655 


5 


325 


30 


min 


California 


2 178 025 


5 


394 


762 


1 010 


5 


341 


86 


min 


Texas 


2 884 762 


7 


194 


984 


1 336 


5 


385 


135 


min 



5 The Route Planning Algorithm - Queries 



Having already computed the boundary graph and 5-reach in the preprocessing 
phase, we now describe a natural simplified two-stage query algorithm based on 
the former. In its cellular stage, as outlined in the introduction, the algorithm 
runs 5-Dijkstra's search until all its branches get saturated at cell boundaries 
(typically, only one or two adjacent cells are searched). Then, in the boundary 
stage, virtually any established route planning algorithm may be used to finish 
the search (cf. Prop. I4.3[) since Bg is a relatively small graph (Table [3]) and is 
free from scope consideration. E.g., we use the standard reach-based A* [T^] . 

Two-stage Query Algorithm (simplified). Let a road network (G,w), a 
proper scope mapping S, an 5-reach r s on G, and the boundary graph 
Bg associated with an edge partition E of G, be given. Assume start and 
target positions s,t £ V(G); then the following algorithm computes, from s 
to t, an optimal si-admissible and saturated walk in G for the scope S. 

1. Opening cellular stage. Let / s C G initially be the subgraph formed by the 
cell (or a union of such) containing the start s. Let r(I s ) C V(I S ) denote 
the actual boundary of I s , i.e. those vertices incident both with edges of I s 
and of the complement (not the same as the union of cell boundaries). 

(a) Run S-reach S-Dijkstra's algorithm (unidir.) on (I s ,w) starting from s. 

(b) Let U be the set of non- s-saturated vertices in r(I s ) accessed in ((Taj. 
As long as U ^ 0, let I s <— I s U Ju where Ju is the union of all cells 
containing some vertex of U, and continue with (|la[). 

An analogical procedure is run concurrently in the reverse network (G R ,w) 
on the target t and It ■ If it happens that I s , It intersect and a termination 
condition of bidirectional Dijkstra is met, then the algorithm stops here. 

2. Boundary stage. Let B s t be the graph created from Bg U {s, t} by adding 
the edges from s to each vertex of r(I s ) and from each one of r(I t ) to t. 
The weights of these new edges equal the distance estimates computed in JTJ. 
(Notice that many of the weights are actually oo — can be ignored — since the 
vertices are inaccessible, e.g., due to scope admissibility or 5-reach.) 

Run the standard reach-based A* algorithm on the weighted graph B s , t (while 
the reach refers back to Bg), to find an optimal path Q from s to t. 

3. Closing cellular stage. The path Q computed in ([2]) is easily "unrolled" into 
an optimal si- admissible saturated walk P from s to t in the network (G, w). 

A simplification in the above algorithm lies in neglecting possible non- 
saturated walks between s,t (cf. Theorem I3.11|) . which may not be found by 
an 5-reach-based search in |T]). This happens only if s, t are very close in a local 
neighborhood wrt. S. 

We provide here a formal description of our query algorithm in pseudocode, 
see Algorithm 15.11 One minor aspect of this algorithm dealing with possible 
"unsaturated" optimal walks in short-distance queries which cannot be properly 
handled by 5-reach-based search, is closely described next. 



We also briefly comment on the use of Dijkstra's and A* algorithms in the 
routing query Algorithm l5.ll While A* is obviously much better in general situa- 
tions, and it is used in the boundary stage, the cellular stage is genuinly different 
in the fact that the search has to get to an s-saturated state in all directions, and 
so A* would be of no help in the cellular stage (the speed-up there is achieved 
by different means). 

One more term related to Definition 14.11 of graph partitioning and cells is 
needed in the algorithm pseudocode: For v £ V(G) let Jg(v) be the set of all 
indices 1 < j < I such that v £ V(Cj). 

Unsaturated Local Search. Relating to line 1 of Algorithm 15. 1[ we now ex- 
plain the necessity of running a separate local search for possible non-saturated 
optimal s—t walks in the road network. As noted in the main paper, it is an 
inevitable limitation of the 5-reach 5-Dijkstra's algorithm that it may not 
search/discover an s-t walk that is optimal but not saturated (cf. Theorem l3.11[) . 
Note that the presence of such a singular walk cannot be simply tested with a 
condition I s n It ^ in the cellular stage of the algorithm, and so a separate 
local search must be used. 

Though, thanks to Def. 13.31 and with well distributed scope mapping, such a 
singular walk may occur only if the distance between s, t is quite short, and then 
even not-so-elaborate route planning algorithms may be used very efficiently. To 
be precise, if a bidirectional search is used, then each direction is searched only 
until s- or i-saturated vertices are reached, which is a quite small subset of the 
whole network. (If these two search spaces do not meet, then no unsaturated s—t 
walk exists.) Furthermore, yet another natural adaptation of the classical reach 
concept — an "anti-5-reach" , can be used to accelerate the search. 

Very briefly, the anti-S -reach of a vertex v is defined as the maximum value 
of minjlP 1 ™!™, [-P^l™} over all pairs of vertices x,y such that there exists a 
non-saturated optimal xy-admissiblc walk P from x to y passing through v, and 
pxv pvy are su b wa l]c S between x,v and between v,y. Again, this parameter 
is easy to compute by brute force in the network. 

5.1 Experimental evaluation. 

Practical performance of the query algorithm has been evaluated by multiple 
simulation runs on Intel Core 2 Duo mobile processor T6570 (2.1 GHz) with 
4GB RAM. Our algorithms are implemented in C and compiled with gcc-4.4.1 
(with -02). We keep track of parameters such as the number of scanned vertices 
and the queue size that influence the amount of working memory needed, and 
are good indicators of suitability of our algorithm for the mobile platforms. 

The collected statistical data in Table |6] — namely the total numbers of 
scanned vertices and the maximal queue size of the search — though, reasonably 
well estimate also the expected runtime and mainly low memory demands of the 
same algorithm on a mobile navigation device. A sample query comparison to 
bidirectional Dijkstra's algorithm is then depicted in Fig. |3l and some statistical 
data are briefly summarized in Table [5j 



Algorithm 5.1 Route Planning Query Algorithm 



Input: A road network (G, w), a proper scope mapping S, the associated <S-reach r s 
(or an upper bound on it), a partition £ of E(G) with the associated boundary 
graph Be (of distance weighting p), and start and target vertices s,t G V(G). 

Output: An optimal si-admissible walk from s to t in this road network (or oo). 
Query(G, to, S, £ , B £ , r s , s, t) 

1: Unsaturated-Local-Search(G, w, S, s, t) //Preliminary "unsaturated" stage. 
Apart from the main query algorithm below, a local search for possible non-saturated 
optimal s-t walks has to be run first (in case of close positions s,t). Details are 
described separately. 

2: I s = Ujej £ ( s ) On I's // Opening cellular stage (from the start). 

3: while I' s / I s do 
4: I' s «- J, 

5: d s ,7r s <- 5-Dijkstra(J s , w rJ - 3 ,<S r / 3 ,s,r' s r /J 
6: 17 {ti G r(/ s ) | VI G Jm(5) \ {oo}. cr £ [u] > uf} 

7: ^<-^.UU jHsm Cj 
8: done 

9: R a <- {v G V(B £ ) | d s [v] < oo} 

10: I t = UjeJ (t) ^* ^ // Opening cellular stage (to the target). 

11: while J t ' / J t do 
12: / t ' <- I t 

13: dt,7Tt <- S-DuKSTRA(7f,w rJt ,<S r / t ,i,r s [G Ji ] r/t ) 
14: U <- {tt G T(/ t ) | W G Jm(S) \ {oo}. <r*[u] > i/f } 

15: Jt«-/«UU jeJe(t0 C, 
16: done 

17: R t <- {v <EV(B £ )\d t [v] <oo} 

18: Bs,t^B £ ; V(B s , t ) <- V(B s ,t) U {s, *} 

19: for all u £ R s do // Boundary stage. 

20: S(S s , t ) <- £(B s ,t) U (s,w); p(s, w) <- d s [u] 

21: done 

22: for all u G Rt do 

23: S(S s , t ) <- E(B s ,t) U (w, t); p(u, t) <- d t [u] 

24: done 

25: d Sjt , («o, «0, • • • ,Ufc, Vfc) «- REACH-BASED-BlDIRECTIONAL-A*(B s ,t, p, S,t) 

Recall that the reach values in this algorithm refer back to Be (and not to B s ,t). 
Therefore, for a vertex v accessed from u, one has to remember which edge e s the 
optimal walk from s to u seen so far starts within B s ,t- The weight w(e s ) must be 
subtracted from d[u] while checking the reach admissibility condition. 

26: for all (m,Vi) G E((uq,vo, . . . ,u k ,v k )) do // Closing cellular stage. 

27: d <s— a cell containing edge (m,Vi) 

28: di,ivi 4- Bidirectional- A*(G°°,w r c~,^i,«0 

29: done 

30: P s , t «- 7r s [tto].7To • • • 7r fc .7rf [v k ] 

31: return d s ,t, P a ,t 



Fig. 3. An example of solving the routing query. The optimal si-admissible walk from 
the start s (green) to the target t (red) is depicted in blue, while the overall optimal 
one in red. The vertices scanned during our query are in light blue, and those scanned 
during bidirectional Dijkstra's search are in light magenta. A brief comparison with 
a map shows that the blue route prefers major highways, making the route more 
comfortable (and perhaps faster in practice). Colorado, CA, United States. 



Table 5. A comparison of our approach with bidirectional Dijkstra's algorithm 
(B-Dijkstra). The US-all (continental part of the US) road network is used. Each row 
lists average for 1 000 uniformly chosen pairs of start and target vertices for distance 
intervals. It could be easily seen that the number of scanned vertices increases only 
a slightly (more vertices of the boundary graph are scanned) in our query algorithm, 
while the number of vertices scanned during bidirectional Dijkstra's search increases 
significantly. 



Distance 


Scanned vertices 


Max. queue size 


(km) 


our approach / B-Dijksta 


our approach / B-Dijksta 


2000 - 3000 


3 225 / 28 139 196 


60 / 4276 


1000 - 2000 


2 854 / 14 821 472 


57 / 2477 


500 - 1000 


2 305 / 2 147 015 


53 / 2011 


0- 500 


1 714 / 692 769 


41 / 1587 



Table 6. Experimental queries for the preprocessed continental US road network (USA- 
all) from Table [3] Each row carries statistical results for 1000 uniformly chosen start- 
target pairs with saturated optimal walks. The table contains average values for the 
numbers of cells hit by the resulting optimal walk, the overall numbers of vertices 
scanned in G during the cellular and in Be during the boundary stage, the maximal 
numbers of elements in the processing queue during the cellular and the boundary 
stages, and the average time spent in the cellular and the boundary stages. We remark 
that this statistics skips the closing cellular stage since the computed walk can be 
"unrolled" inside each cell on-the-fly while displaying details of the route. 



Query distance 


Hit 


Scanned vertices 


Max. queue size 


Query time 


(km) 


cells 


cellular / boundary 


cellular / boundary 


(ms) 


3000 - oo 


277 


1 392 / 3 490 


60 / 58 


8.2 + 29.8 


2000 - 3000 


139 


1 411 / 1 814 


64 / 52 


7.9 + 26.9 


1000 - 2000 


57 


1 343 / 1 511 


57 / 49 


7.7 + 22.8 


500 - 1000 


25 


1 113 / 1 192 


53 / 38 


8.1 + 19.0 


0- 500 


10 


998 / 716 


41 / 34 


6.9 + 16.1 



6 Conclusions 

We have introduced a new concept of scope in the static route planning problem, 
aiming at a proper formalization of vague "route comfort" based on anticipated 
additional metadata of the road network. At the same time we have shown 
how the scope concept nicely interoperates with other established tools in route 
planning; such as with vertex-separator partitioning and with the reach con- 
cept. Moreover, our approach allows also a smooth incorporation of local route 
restrictions and traffic regulations modeled by so-called maneuvers [17] . 

On the top of formalizing desired "comfortable routes" , the proper mixture of 
the aforementioned classical concepts with scope brings more added values; very 
small size of auxiliary metadata from preprocessing (Table [3]) and practically 
very efficient optimal routing query algorithm (Table [6]). The small price to 
be paid for this route comfort, fast planning, and small size of auxiliary data 
altogether, is a marginal increase in the weight of an optimal scope admissible 
walk as compared to the overall optimal one (scope admissible walks form a 
proper subset of all walks) . Simulations with the very basic scope mapping from 
Table Q] reveal an average increase of less than 3% for short queries up to 500km, 
and 1.5% for queries above 3000km. With better quality road network metadata 
and a more realistic utility weight function (such as travel time) these would 
presumably be even smaller numbers. 

At last we very briefly outline two directions for further research on the topic. 

i. With finer-resolution road metadata, it could be useful to add a few more 
scope levels and introduce another query stage (s) "in the middle". 

ii. The next natural step of our research is to incorporate dynamic road network 
changes (such as live traffic info) into our approach — more specifically into 
the definition of scope-admissible walks; e.g., by locally re-allowing roads of 
low scope level nearby such disturbances. 
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